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Abstract 



' We consider the Einstein- Yang- Mills Lagrangian in a (4 -|- n)-dimensional space-time. Assuming 

■ the matter and metric fields to be independent of the n extra coordinates, a spherical symmetric 

i Ansatz for the fields leads to a set of coupled ordinary differential equations. We find that for 

I n > 1 only solutions with either one non-zero Higgs field or with all Higgs fields constant exist. We 

construct the analytic solutions which fulfill this conditions for arbitrary n, namely the Einstein- 
Oh. Maxwell-dilaton solutions. 



X 



We also present generic solutions of the effective 4-dimensional Einstein- Yang-Mills-Higgs-dilaton 
model, which possesses n Higgs triplets coupled in a specific way to n independent dilaton fields. 



. These solutions are the abelian Einstein-Maxwell- dilaton solutions and analytic non-abelian solu- 

tions, which have diverging Higgs fields. In addition, we construct numerically asymptotically flat 
and finite energy solutions for n = 2. 
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I. INTRODUCTION 



In an attempt to unify electrodynamics and general relativity, Kaluza introduced an 
extra - a fifth dimension l| and assumed all fields to be independent of the extra dimension. 

n 

Klein j2| followed this idea, however, he assumed the fifth dimension to be compactified 
on a circle of Planck length. The resulting theory describes 4-dimensional Einstein gravity 
plus Maxwell's equations. One of the new fields appearing in this model is the dilaton, a 
scalar companion of the metric tensor. In an analogue way, this field arises in the low energy 
effective action of superstring theories and is associated with the classical scale invariance 
of these models 0. 

When studying spherically symmetric solutions in higher dimensional systems, two pos- 
sible approaches seem possible: a) to assume the solutions to be spherically symmetric in 
the full d dimensions or b) to assume the solutions to be spherically symmetric only in 4 
dimensions. The solutions obtained in the first approach are surely important at very high 
energies, i.e. very early stages of the universe. The second approach, which assumes the 
extra dimensions to be "spectator" is of importance for physics in the universe today. 

Volkov recently followed the idea of extra dimensions and constructed a (4+1) dimensional 
Einstein- Yang-Mills (EYM) system j^. It turned out that the EYM particles are completely 
destroyed by gravity, but that, however, particle-like solutions (so-called "EYM vortices") 
exist if one assumes the fields to be independent of the extra coordinate. The system then 
reduces to an effective 4-dimensional Einstein- Yang-Mills-Higgs-dilaton (EYMHD) model, 
which was studied in detail in P| . Through the dimensional reduction, one Higgs triplet and 
a dilaton field appears in the model. 

n 

In this paper, we study the system of 4] in (4 + n) dimensions, assuming all fields to 
be independent of the n extra coordinates. We give the model and equations of motion 
in Section II. The solutions are spherically symmetric in 4 dimensions, while the extra 
dimensions (which are "spectator") are associated with a Ricci fiat manifold. It turns 
out that in n > 1 dimensions only solutions with one non-zero Higgs field or with all 
Higgs fields constant exist. We give the analytic solutions which fulfill this condition for 
generic n, namely the abelian Einstein-Maxwell-dilaton (EMD). In Section HI, we present 
the 4-dimensional effective Einstein- Yang-Mills-Higgs-dilaton (EYMHD) model, which has 
n independent Higgs and dilaton fields. Since in this model, no constraint on the Higgs 
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fields arises, we construct generic solutions. We give the analytic solutions available, namely 
the abelian Einstein-Maxwell-dilaton (EMD) solutions and "non-abelian" solutions with 
diverging Higgs fields. Finally, we present our numerical solutions of the effective model for 
n = 2. We give our conclusions in Section IV. 

II. THE (4 + n)-DIMENSIONAL EINSTEIN- YANG-MILLS MODEL 

The Einstein- Yang-Mills Lagrangian in d = {4 + n) dimensions is given by: 

with the SU(2) Yang-Mills field strengths F^^^ = OmA^ - OnAI^ + eahcA\jA% , the gauge 
index a = 1,2,3 and the space-time index M = 0, (4 + ra) — 1. G'(4+„) and e denote 
respectively the (4 + n) -dimensional Newton's constant and the coupling constant of the 
gauge field theory. G'(4+„) is related to the Planck mass Mpi by G(4+n) = Mp^*-^'''"^ and 
has the dimension of [length]". In the following, we denote the coordinates X(3+fc) by t/k with 
= 1, n. 

If both the matter functions and the metric functions are independent on yk, the fields 
can be parametrized as follows: 

n 

g^^^dx^dx"" = e~-gl^Jdx>'dx'' + e^^" {dy^f , /i, z/ = 0, 1, 2, 3 (2) 

k=l 

with 

n 

s = Ec. (3) 

k=l 

and 

n 

A%dx^^ = A^^dx" + J2 ^tdy^ ■ (4) 

k=l 

g^^ is the 4-dimensional metric tensor and the Q and j = l,...,n, play the role of 
dilatons and Higgs fields, respectively. 

Note that in the case n = 1, the above parametrization coincides with the one in jj]. 
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A. Spherically symmetric Ansatz 



For the metric the spherically symmetric Ansatz in Schwarzschild-like coordinates reads 
ds^ = glf^dx^'dx" = -A\r)N{r)dt^ + N-\r)dr^ + r^dO^ + sin^ ed^^ , (5) 



with 



iV(r) = l-^. (6) 



In these coordinates, m{oo) denotes the (dimensionful) mass of the field configuration. 
For the gauge and Higgs fields, we use the purely magnetic hedgehog ansatz 0| : 

A/ = A" = , (7) 

Ae'' = {l-K{r))e^\ A^'' = -{1 - K{r)) sin Oee", (8) 

^'^ = c,vH,{r)er'' , j = l,...,n, (9) 

where w is a mass scale, while Cj are dimensionless constants determining the vacuum ex- 
pectation values of the Higgs fields = CjV. In absence of a Higgs potential these have 
to be set by hand. Finally, the dilatons are scalar fields depending only on r : 

Q = Cj{r) , k = l,...,n. (10) 

With these conventions, the non- vanishing components of the (4 + ra)-dimensional energy 
momentum tensor read: 



^0 = -e- [T.yk + Ck) + B + Dj , 
Tl =-e^(^J2i-^k + Ck)-B + D^ , 
Ti =Ti = -e^(pA,-D^ , 

T/ =-e-(^-2{A,+C,) + j2iAk + Ck) + B + Dj , j = l,...,n, 

/ N \ 
= e^e-(^^+«'=) [i^H'^K + -^H.H,^ . J^k , (11) 



where we use the abbreviations 



A, = \e-^^^N{H'^f , C, = e-'^^K^H]^^ (12) 
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and 



B-e^^^N(Kr , D-e^^^(K^-ir (13) 



The radial variable x = evr was introduced and the prime denotes the derivative with respect 
to X. 



B. Equations of motion 

Defining the coupling constant a — v^Jc'i^^) and the mass function ji — evm we obtain 
two first order equations for the two metric functions 

ix' = a^x'' (fli^k + Cfc) + S + L>j + ^TVx^e , (14) 

A' = a^AxiY.e-^'^'iH'^f + 2e~\ \ + xAQ (15) 
\fe=i ^ / 



with 



E(0' + ?Ea^'| • (16) 



4 \ 3 

^ \it=l k>k' 



The components of the Einstein equations related to the extra dimensions lead (after suitable 
linear combinations) to equations for the n dilaton fields : 

d-l 



{x'ANQ' = a'Ax'e 



2 

o 

M=0 



, j^l,...,n. (17) 



The field equations for the gauge field and the n Higgs fields read respectively 

(e^ANKj ^A(^e^K{K'' -l) + f2 e-^^'K Hi] , (18) 

k=i J 

(e-'^^^x'^ANH'j)' ^2Ae-^'^^K^Hj , j^l,...,n. (19) 

Finally, since the off-diagonal components of the Einstein tenor vanish, we obtain an extra 
constraint on the fields from the j A;- components of the energy-momentum tensor: 



(20) 



Clearly, these equations are symmetric under the simultaneous exchange Q'lHji Cj^Hj. 
Further, they are invariant under the rescaling 

Hj X-^Hj , X ^ Xx , X^a , X e IR . (21) 
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The constraint ()20j) will only be fulfilled for specific cases. With our parametrisation, 
there are two possibilities: (i) only one non-zero Higgs field and arbitrary gauge field K{x) 
or (ii) constant Higgs fields and K{x) = 0. 

Since we are interested in generic solutions, we go back to the 4-dimensional effective ac- 
tion motivated by this model in which, of course, off-diagonal terms of the energy-momentum 
tensor don't appear. Before we discuss this effective model, we give the solutions with con- 
stant Higgs fields available in the "full" model, namely the Einstein-Maxwell-dilaton (EMD) 
solutions. 

Let us remark that there are several possibilities of modification of the "full" model 
considered here, which could eventually lead to solutions with generic Higgs and gauge 
fields, namely we could 

(i) add a non-diagonal term of the form f^Qdx^dx^ in the metric, this leads to one extra 
equation for the /sg function, but it has to be checked that this function stays regular, 

(ii) start with an SU(N) gauge group and choose the different Higgs fields in orthogonal 
SU(2) subalgebras of the Lie algebra of SU(N). 

C. Einstein-Maxwell-Dilaton (EMD) solutions 

Assuming the matter fields to be constant: 



we find that the above equations ((13)), (fTKjl . (|T7)l admit exact solutions which are related to 
the Einstein-Maxwell-Dilaton (EMD) solutions P]. Here, we will only discuss the extremal 
case, which corresponds to a solution with horizon at the origin. 
The extremal EMD solutions have unit magnetic charge and mass 



K{x) = , Hj{x) = Cj , j = 1, ...,n 



(22) 



and 



0(x) = C(x) Vj, 



(23) 




(24) 



The value of the metric component N{x) at the origin x = reads: 




(25) 
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and the dilaton field is: 



CEMD{x) = ^\n[l-^) (26) 



with 



III. THE EFFECTIVE 4-DIMENSIONAL EINSTEIN- YANG-MILLS-HIGGS- 
DILATON (EYMHD) MODEL 

As in the 5-dimensional case P] (in our notation n = 1) the equations given in the 
previous section (apart from the constraint ((201)) can equally well be derived from an effective 
4-dimensional Einstein- Yang-Mills-Higgs-Dilaton (EYMHD) Lagrangian. The Lagrangian 
density for the matter fields then reads: 

1 " 1 

^ k=i 2 

J2 d^^kd^^k + I E df.'^kd^'^k' I (28) 

^fc=l k>k' 

with 

r = E^. . (29) 

k=l 

The kinetic part in the dilaton fields could be diagonalized, however, we find it more conve- 
nient to leave it in the form above which reveals the symmetry $j 

The Lagrangian (PHj) is then coupled minimally to Einstein gravity according to the full 
action 

S = Sq + Sm 

= J 4-9^^ (^G + Lm) d^x (30) 

where Lq = R/ (167rG'4), R is the Ricci scalar and G4, is the 4-dimensional Newton's constant. 

Note that the dilaton fields are coupled by an independent coupling constant n to the 
gauge and Higgs fields. In this respect, the dilatons here are treated as independent scalar 
fields, while in the action they appear as parts of the metric tensor. 

After the rescaling 

1 I — 
'^j = vibj , K, = — , a = v\ Gi (31) 

V ^ 



the resulting set of equations only depends on the coupling constants a and 7. We refrain 
from giving the explicit form of the equations here, but refer the reader to 5| for the case 
n = 1. 

Note that the equations (HH), (fT^ . (fTTj) . (HHj), (fT^ become equivalent to the field equa- 
tions associated to (|3(jp by using the same Ansatze for the 4-dimensional metric, the gauge 
and Higgs fields, but by identifying 

Cj = = 2k^j , = 87^ . (32) 

Remarkably, this identification turns out to be independent on n. 

This model now has solutions for which all n Higgs fields can be non-constant. First, 
we will present the 4-dimensional effective counterparts of the abelian Einstein-Maxwell- 
dilaton solutions. Then, we will present analytic non-abelian solutions with non-constant 
but diverging Higgs fields, which are not available in the full model since they don't fulfill 
the constraint (jSOl)- Finally, we will give the asymptotically flat, finite energy solutions for 
n = 2, which we construct numerically. 



A. Einstein-Maxwell-Dilaton (EMD) solutions 

Assuming the matter fields to be constant: 

K{x) = , Hk{x) = Cfc , k = l,...,n (33) 

and 

Mx)=^ix) Vfc, (34) 

we find that the equations associated to the effective action (PUI) admit exact solutions which 
are related to the Einstein-Maxwell-Dilaton (EMD) solutions [3]. As in the previous section, 
we will only discuss the extremal case. Note that now, of course, the solution depends on 
both a and 7. 

The extremal EMD have unit magnetic charge and mass 

A^oo 1 



where 



(35) 



■? - J^f ■ (36) 
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The value of the metric component N{x) at the origin x = reads: 



and the dilaton field is: 



NEMniO)=(-^J (37) 



^emd{x) = ( — — ) In fl - ^ ) (38) 



.2 + nJ a2 + ^2-^- X 

with 



B. Analytic, non-abelian solutions 



Very similar to the n = 1 case |4| non-abelian, analytic solutions exist in this system. 
If we assume the limit = 87^ (which, of course, is the most interesting one), we find 
solutions of the form: 

N{x) = ci , A = x(^+"^) , Mx) =ip2ix) = ... ^il^ix) = ^L\n(— 

/ „ \ 2/n 

H^{x) = H2{x) = ... = H{x) = ,/FJ-) , K{x) = ^ (40) 



-C3/ 

where Ci, C2, C3 and C4 are constants depending only on n, g and a and are given by: 

-2q^rr'in + l) + qn^i2n^ + bn + &)-n'^ 2, , , 

^1 = ^ o / ^^^ \A(o ^ ' C2 = - l + 2a2c4 , 

8gri(n + 1) + 4(2 — n) n 

2 -4g^a2(n + 1) + 2q^a^{2n'^ + 7n + 8) + Aqa'^{-n^ - 3n - 3) + 2a^n 

~ 2g2n(-n2 - 3n - 2) + q{2n^ + Sn^ + Gn^ + 2n - 4) - ' 

C4 = ^ , (41) 
and g satisfies the equation 

g^(2ra^+4n^+2n2)+g2(-2r2^-3n^+ra^+10ra2+8r2)+g(3n^+2n^+3ra2+6r2+8)-ra3+n2-2n = . 

(42) 

Solving this condition for q numerically for various values of n, we find one real positive 
solution for n = 1, 2, 3 and three positive solutions for n = 4, 5, 6, 7. However, when more 
than one solution for q are available, it turns out that only one of the solutions is acceptable 
since for the others some of the remaining parameters (e.g. C4 or C3) become negative. In 
Table 1, we give the acceptable solutions and their numerical values for the parameters Ci, 
C2, C3, C4: 



Table 1: Values of q and constants ci, C2, C3, C4 for different n 



n 


= 


Cl 


C2 




a^C4 


1 


0.293 


0.016 


5.246 


1.126 


0.811 


2 


0.294 


0.064 


1.702 


1.301 


0.351 


3 


0.280 


0.102 


0.982 


1.297 


0.237 


4 


no real solution 










5 


0.249 


0.149 


0.495 


1.581 


0.118 


6 


0.236 


0.151 


0.554 


0.476 


0.331 


7 


0.225 


0.086 


2.402 


0.036 


3.704 



C. Numerical results for n = 2 

The 4-dimensional effective theory for n = 1 was studied in j^. 

Here, we will put the emphasis on the case n = 2 and again discuss the solutions in the 
limit = 37^. In this model, two Higgs fields naturally occur. Theories involving two 
Higgs fields are also interesting from the viewpoint of the supersymmetric extension of the 
Standard model. 

The case n = 2 involves two parameters only, namely a and p = C2/C1. As discussed 
previously in this paper, it is sufficient to consider p > 1. 
The solutions can then be characterized by their mass : 

1 /i(c)o) 



M 



c\ + cl « 



2 ' 



(43) 



where the first factor is extracted in such a way that the mass of the solution in the a = 
limit is normalized to 1, corresponding to the BPS monopole. 

The values '?/^i,2(0) will appear to be useful to further characterize the solutions. The 
deviation of the solutions from the flat space solutions can be "measured" by the values of 
A{x) at the origin, A(0), and by the minimum of N{x), say Nm- Since A{x) always increases 
monotonically, A{0) also represents the minimum of the function. 

In the following, we will discuss our numerical results for two different and fixed values 
of p and varying a. 

We again discuss the solutions in the limit of = 37^. 
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1. Boundary conditions 



We will study globally regular, asymptotically flat solutions of the system above. This 
implies the following boundary conditions: 

ir(0) = l, i/,(0)=0, a,.^,U=o = 0, M0) = (44) 

at the origin and 

K{oo) = , Hj{oo) = Cj , ^j(oo) = , A{oo) = 1 , (45) 

at infinity. We assume Cj 7^ in order for the solutions to have a magnetic charge. Using 
the symmetries discussed above, we see that we can set min{cj,j = l,...,n} = ci = 1 
without loosing generality and study the equations with varying a and 02,03, ...,c„ > 1. 
Of course, in the case ci = C2 = ... = = 1 the equations are symmetric and we expect 

Hi{x) = H2{x) = ...Hn{x) = H{x) , ^i(x) = ij2{x) = ... = i^nix) = ^'(x). 



2. p=l 

In this case, the two Higgs functions Hi{x) and H2{x) are equal. This, of course, implies 
immediately that also ipi{x) = 'ip2{x)- 



In the fiat limit {a = 0) the solution is the BPS monopole 10|,lll| with A{x) = N{x) = 1, 
ipiix) = ^'2(2;) = and the matter functions K{x), Hi{x) = H2{x) have the well known 
BPS profiles. 

Increasing a, our numerical results reveal that the solution gets progressively deformed 
by gravity. The solutions form a branch on which the mass M dimnishes with increasing a. 
The same holds true for the values ^4(0), Nm- 

The function N{x) indeed develops a minimum which becomes deeper while gravity 
increases. At the same time the dilaton functions ipi.2{x) are non trivial for a ^ 0. The 
value of V'i,2(a; = 0) is negative and decreases with increasing a. 

Our results are illustrated in Fig. ^ The figure also demonstrates that the branch of 
gravitating solutions does not exist for arbitrarily large values of the parameter a. Indeed, 
we find that solutions cease to exist for a > amax '^mLx ^^us constitutes a maximal value 
of a. We find that this critical value depends on the parameter p and is given here by: 

a(^L ^ 0.75 for p = 1 . (46) 
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Our numerical analysis further reveals that for a < a^^^ another branch of solution 
exists. The mass of the solutions on this second branch is higher than the mass of the 
corresponding solutions on the first branch. This is demonstrated in Fig. |21 

Our numerical results further suggest as illustrated by the figures that the solutions on 
the second branch stop to exist at some local minimal value a = amm- 

To be more precise, the further evolution of the branches is involved : indeed several small 
branches exist for a e a^;^ax^] with a smaller and smaller extend in a. Nevertheless, 

it appears clearly from the figures that for a acr, where acr '■= «min = '^mtx^ some 
k, the value ^4(0) tends to zero, while the values |V'i,2(0)| increase considerably and likely 
become infinite. Note, however, that the additional branches are not visible in the mass plot 
(see Fig|21), because the numerical values are very close to those of the second branch. The 
existence of several branches was also noticed for the n = 1 system in iM and in another 

n 

Einstein- Yang-Mills model in 5 dimensions |12l |. 
We find numericalUy: 

aSL~ 0.240 , a^L^ 0.294 , a,, ^ 0.285 (47) 

At the same time, the value Nm stays strictly positive and tends to the value given by 
()40|) in the critical limit. We find numerically that ~ 0.065 for a ~ acr, while Ci ~ 0.064 
for n = 2. The function N{x) reaches this minimum at a value of x, which tends to zero for 
a acr- At the same time, the function N{x) becomes nearly fiat on a plateau surrounding 
this minimum. 

The function K{x) starts to develop oscillations around the value ^/g = 0.294. Processing 
on the branches, the number of oscillations increases and becomes infinite in the critical limit. 



Very similar to what was observed in 



the Yang-Mills domain walls get disconnected 



from the outside world by an infinitely long throat in the strong gravity limit. 



3. p = 2 

In order to understand the influence of two Higgs fields with different expectation values, 
i.e. p 7^ 1, on the domain of existence of the solutions, we studied in detail the case p = 2. 

Varying a, the existence of several branches of solutions is qualitatively very similar to 
the case p = 1. The numerical analysis, however, reveals that for p > 1, the solutions exist 
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on a smaller interval of the coupling constant a in comparison to the case p = 1. Namely, 
we find: 

~ 0.496 , aSn ~ 0.150 , a^^L = 0-185 , ^ 0.179 for p = 2 . (48) 

In the critical limit for a ~ Ocr, we find 

M ^ 2.05 , Nm ~ 0.068 (49) 

One of the apparent differences is that the two dilaton fields are non-equal. The values of 
?/'i(0) and ^"2(0) are superposed on Fig. El The difference ^"2(0) — '?/'i(0) is also shown and 
indicates that in the critical limit, the two dilaton profiles remain significantly different close 
to the origin. 

The profiles of the function K{x), N{x), ipii^) and ip2{x) are shown in Fig.EJfor a = 0.179. 
This solution is already close to the critical solution. Note that the function K{x) starts to 
develop oscillations completely analogue to the case p = 1. 

The numerical results suggest that in the critical limit a solution of the analytic, non- 
abelian type described in the previous subsection is reached. The question how this solution 
looks like for p 7^ 1 will be addressed in a future publication. 



IV. CONCLUSIONS 

Both string theories [l^ as well as so-called "brane worlds" jl^ , which assume the Stan- 
dard model fields to be confined on a 3-brane (that is embedded in a higher dimensional 
space-time) have enhanced the idea that space-time possesses more than four dimensions. In 
the former, the extra dimensions are compactified on a scale of the Planck length, while in 
the latter they are non-compact. Non-perturbative, classical solutions of field theory models 
certainly play a major role in these theories. 

It is therefore natural to investigate the classical solutions of higher dimensional Einstein- 
Yang-Mills theory. In this paper, we have investigated Einstein- Yang-Mills theories in 
1+3+n dimensions, in which the matter and metric fields are chosen to be independent 
of the extra n coordinates. Dimensional reduction then leads to n dilaton fields coupled in- 
dividually to n Higgs fields. These n Higgs fields can have independent vacuum expectation 
values. 
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The corresponding spherically symmetric equations admit several types of solutions: (a) 
vacuum solutions, (b) embedded abelian Einstein-Maxwell-dilaton solutions, (c) non-abelian 
solutions with diverging Higgs fields and (d) fully non-abelian magnetically charged solu- 
tions. The vacuum solutions (a) are trivial, while (b) and (c) were constructed analytically 
for the case when all Higgs fields have identical expectation values. The solutions of type 
(d) had to be constructed numerically. 

We believe that type (b) and (c) solutions can also be constructed analytically for generic 
values of the Higgs's expectation values and we plan to reconsider this problem in a future 
publication. 

It would also be interesting to study this model for larger gauge groups or to adopt 
a different compactification scheme of the codimension space, which here was chosen to be 
X X . . . . Following the investigation for n = 1 Il5|l, the model including a cosmological 

n 

constant is presently under investigation. In [15|] the introduction of a cosmological constant 
in the (4+l)-dimensional model led to a Liouville type potential in the effective 4-dimensional 
theory. 
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FIG. 1: The values ^(0), Nm-, —"0(0) ^-re shown as functions of a for p = 1. The indices "1", "2", 
"3", "4", respectively correspond to the 1., 2., 3. and 4. branch of solutions (see also Fig. I^J. 



16 




FIG. 2: The mass M of the solutions is shown as function of a for p = 1 and p = 2. "1" , "2" and 
"3", respectively denote the 1., 2. and 3. branch of solutions. 
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FIG. 3: The values of — ■i/'i(0) and — 1/'2(0) as well as the difference ■i/'2(0) — ipi{0) are shown as 
functions of a for p = 2. "1", "2", "3" and "4", respectively denote the 1., 2., 3. and 4. branch of 
solutions. 
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FIG. 4: The profiles of the functions iIji{x), 'tp2{x), K{x), A{x) and N{x) are shown for p = 2 and 
a close to the critical value CKcr ~ 0.179. 
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